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Abstract. We study the asymptotic behavior of the partition function and the correlation 
kernel in random matrix ensembles of the form | det —tl )where M 
is an n X n Hermitian matrix, a > —1/2 and t G M, in double scaling limits where n —l oo 
and simultaneously t —1 0. If t is proportional to Ijn^, a transition takes place which 
can be described in terms of a family of solutions to the Painleve V equation. These 
Painleve solutions are in general transcendental functions, but for certain values of a, they 
are algebraic, which leads to explicit asymptotics of the partition function and the correlation 
kernel. 
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1 Introduction 


We consider a random matrix model on n x n Hermitian matrices, defined by a probability 
measnre of the form 

^ |det (M^for a >-1/2, t G M. (1.1) 

Here dM = n dMu n dReMjjdlmMjj is the Lebesgue measure on n x n Hermitian 

l<i<n 

matrices and = Zn{t,a, V) is a normalizing constant. The confining potential H: M —>■ M is 
a real analytic function which increases sufficiently fast asymptotically, namely 


V{x) 


lim ^— / 2 I 
log(x^ -k 1) 


( 1 . 2 ) 


The measure (1.1) is invariant under unitary conjugation and induces a probability measure 
on the eigenvalues Xi, i = 1,2,... ,n, given by 


^ n {Xj - Xif'^Wn{Xj)dXj, Wn{x) = -t\°' 

l<i<j<n j=l 


(1.3) 


*This paper is a contribution to the Special Issue on Asymptotics and Universality in Random Matrices, 
Random Growth Processes, Integrable Systems and Statistical Physics in honor of Percy Deift and Craig Tracy. 
The full collection is available at http://www.emis.de/journals/SIGMA/Deift-Tracy.html 
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where 

n 

{Xj - Xif'^Wn{Xj)dXj (1.4) 

l< 2 <j<n j=l 

is the partition function associated to the ensemble. It can alternatively be written as the Hankel 
determinant 

Zn{t,a,V) = nldet ( f x'^~^^Wn{x)dx'\ 

\J—oo J i1 

While a is kept fixed and does not depend on the matrix size n, we will be particularly 
interested in double scaling limits where t is n-dependent and approaches 0 as n —)• oo. In such 
a case, the weight function Wn has the feature that it has two merging singularities approaching 0 
as n —)• oo. If t > 0, the two singularities ±\/t are real; if t < 0, the singularities are 

purely imaginary. Random matrix ensembles with such merging root-type singularities have 
applications in fractional Brownian motion with Hurst index H = d [19], and in quantum 
chromodynamics [1, 6]. The partition function Zn appears in the study of impenetrable bosons 
in a harmonic well as the one body density matrix (up to a constant) [17]. 

Although the singularities have no effect on the macroscopic limiting density of the eigenva¬ 
lues as n ^ oo, they do have an effect on the partition function Z^ and on the microscopic corre¬ 
lations between eigenvalues near the origin. We will obtain asymptotics for the partition function 
and for the eigenvalue correlation kernel near the origin, in double scaling limits where n —)• oo 
and t —7- 0 simultaneously. The asymptotics for the partition function are described in terms of 
solutions to the fifth Painleve equation; the limit of the correlation kernel can also be expressed 
in terms of functions related to the Painleve V equation. As t —)• 0 at a fast rate, the limiting 
kernel degenerates to a Bessel kernel, and as t —)• 0 at a slow rate, we recover the sine kernel. 


(1.5) 


Zin — 


Zn{t,a, V) = 


/ oo pc 

-oo J —C 


1.1 Statement of results 

The macroscopic large n behavior of the eigenvalues is described by the limit of the mean 
eigenvalue distribution as n —)• oo, which we denote by /ry, and which is independent of t 
and a, for t, a bounded. It is characterized [7, 10] as the unique equilibrium measure /iy which 
minimizes 

log^—— -d^i{x)d^i{y) + I V{x)dy,{x), (1.6) 

\x-y\ J 

over the space of Borel probability measures y on M. The equilibrium measure is alternatively 
characterized by the Euler-Lagrange variational conditions [30] 

2 J log \x — y\dyviy) = ^{x) — i, for x £ supp yv: 

2 J log \x — y\dyv{y) < V(x) — £, for x G M\ supp yv, 

for some constant £ depending on V. 

Throughout the paper we will assume that V is such that the equilibrium measure is supported 
on a single interval [a, b]. Then, the equilibrium measure can be written in the form [10] 

dyv{x) = 'ipv{x)dx = \/ {b — x){x — a)h{x)dx, x£[a,b], (1-9) 

with h a real analytic function. We require V to be such that the following generic conditions 
hold: the density V^y is strictly positive on (a, 6), it vanishes like a square root at the end¬ 
points a, b (in other words, h{x) > 0 for x G [a,b]), and the variational inequality (1.8) is strict 


(1.7) 

( 1 . 8 ) 
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on M\[o, h]. If all those assumptions hold, we say that V is one-cut regular (see [26] for a general 
classification of singularities). In addition, we require zero to be contained in the interior of the 
support of /xy, i.e., a < 0 < 6. 

If zero were outside of the support of fiy, no eigenvalues are expected near the origin for 
large n, and the singularities of the weight will not play a significant role for large n. The case 
where 0 is an edge point of the support of V is also interesting, but will not be studied here - we 
refer to [22] for results about the local eigenvalue correlations when a singularity lies near a soft 
edge, and to [33] when a singularity approaches a hard edge in a modified Jacobi ensemble. 


Painleve V functions 

In order to describe the asymptotic behavior in the matrix model (1.1) for n large and t small, 
we need a special 1-parameter form of the more general 4-parameter Jimbo-Miwa-Okamoto 
(T-form of the fifth Painleve equation (see [23, 24], and [18] for the explicit relation to the usual 
Painleve V equation), which is given by 

(scj''(s))^ = (<t(s) — sa'{s) + 2{a{s))‘^ + 2aa\s))'^ — 4:{a'{a'(s) -t- a)^. (1-10) 

For a > —1/2, we are interested in two particular solutions ct+ and a~ to this equation. The 
first one, cr~(s), is analytic and real for s G (0,-t-oo), and it has the asymptotics 




-I- o(l). 


s —)• 0, 
s —)> -|-oo. 


( 1 . 11 ) 


where T is Euler’s T-function. The existence of such a solution was proved in [4]. The second 
solution (T^(s) is analytic and such that o'q(s) ^ is real-valued for s G —iM>o, satisfying in 
addition: 


{ a^-|-o(l), s —)> — iO"'', 

o? as 

Y - Y + ^(l^l )> s^-ioo. 

Its existence was proved in [5] . It should be noted that uniqueness of the solutions satisfying 
the above asymptotic conditions was not proven in [4, 5]. Those Painleve solutions a~ and 
can be characterized in terms of Riemann-Hilbert (RH) problems; we will do this in Section 2 
for a~ and in Section 3 for (t+. For general values of a > —1/2, are transcendental functions, 
but there are special values of a where they are algebraic. For (T“) is the case if a is an 
integer, while for only if a is an even integer. For those values of a, can be constructed 
via a recursive procedure, involving Schlesinger type transformations [14], which we will explain 
in detail in Sections 2.2 and 3.2. 

For a = 0,1, 2, the expressions are simple and we have 


Cj^(s) = 0, 

(1.12) 

("l-e.-l' 

(1.13) 

^ , —2 -|- e*(2 — 2s -|- s^) 

^2 IS) - s—— eT(2 ■'+ s^y + ' 

(1.14) 


Asymptotics for the partition function 

Theorem 1.1. Let V be real analytic, satisfying (1.2) and one-cut regular with 0 in the interior 
of the support [a,b] of the equilibrium measure pLy Denote zq = y/t when t is positive and 
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zo = iy/^ when t is negative. Define 

/ zo 

h{s) {{s — a){b — ds, 

-20 

with h defined in (1.9) and where the square root is taken to he positive at 0 and analytic in 
a neighbourhood of 0, in such a way that Sn,t is positive when t < 0 and negative imaginary 
when i > 0. We assume that a > —1/2 and t < 0 or that a > 0 and t > 0. Then, the partition 
function Zn{t) defined in (1.4) satisfies 

__ _ + 

log Zn{t, a, V) = log Zn{0, a, V) + - - ds -\ -^ 

Jo s 2 

+ "^(Vizo) + V{-zo) - 2F(0)) + 0{\t\^/^) + O(n-i), (1.15) 

with uniform error terms as n —?• oo and t —)• 0, > 0. 


Remark 1.2. If we set 'Sn,t = —4:7rinzofiv{0), we have finy = Sn,t (1 + 0(|^|)) as n —>■ oo, t —>■ 0. 
If t > 0, it follows from equation (1.36) below that one can replace Sn,t by 'Sn,t in (1.15) without 
modifying the error terms. Doing the same when t < 0 induces an error term of order 0(n|t|^/^). 

Remark 1.3. We believe that the expansion (1.15), with error term o(l), holds for t > 0 and 
—1/2 < a < 0 as well, but because of technical obstacles, we do not prove the result in this case. 

Remark 1.4. The asymptotic behavior for Zn{t,a,V) is described in terms of the Painleve V 
solutions and in terms of the partition function Znfi),a,V), which corresponds to a weight 
function which has only one singularity at 0. Asymptotics for Zn{t),Oi, V) were obtained in [25] 
in the case where V is quadratic, an important special case which we will return to later in this 
section. For general V, we have not been able to find explicit asymptotic expansions in the 
literature for logZ„(0, a, V), up to decaying terms as n —>■ oo. 


Remark 1.5. By (1.12)-(1.14), the integral appearing in (1.15) can be computed explicitly for 
a = 0,1, 2. We have 





ds = 0, 


(1.16) 


J (^4_(^- Ids = log ^2 sink ^-log Sn,t, (1.17) 

J -ds = log (^4sinh^ ^ - 41ogSn,t+ logl2. (1.18) 

Remark 1.6. Our proofs of Theorems 1.1 and 1.7 are based on an asymptotic analysis of the 
orthogonal polynomials on the real line defined by (1.19), using RH methods [7, 11, 12], The 
asymptotic analysis for those orthogonal polynomials shows similarities to the one performed 
in [4, 5] (following the general method from [8]) for orthogonal polynomials on the unit circle 
with respect to weights with merging or emerging singularities, where they were used to obtain 
asymptotics for Toeplitz determinants. Nevertheless, there is no known way to directly derive 
the asymptotics for the Hankel determinant (1.5), with respect to the n-dependent weight Wn 
supported on the full real line, from the asymptotics for Toeplitz determinants. 


Asymptotics for the correlation kernel near 0 

Denote by p^^'^ the degree j orthonormal polynomials with respect to the weight Wn{x) (defined 
in (1.3)), characterized by 



Wn{x)dx 


djk — 


0 for j / k, 
1 for j = k, 


(1.19) 
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with positive leading coefficient The eigenvalue density (1.3) characterizes a determinantal 
point process with correlation kernel given by 

n—1 

Kn{x,y) = ^/wn{x)Wn{y)'^vf\x)vf\y), 

j=0 


or alternatively 


Kn{x,y) = sj Wn{x)'Wn{y) 


4-1 Pn^ {x)Pn\{y) - Pn-l{x)Pn^ (?/) 


(n) 

Kn 


x-y 


( 1 . 20 ) 


by the Christoffel-Darboux formula. Note that the orthogonal polynomials and the correlation 
kernel Kn depend on the parameter t in the weight Wn- This paper is concerned with the 
microscopic large n behavior of the eigenvalues near 0. We will study the scaled correlation 
kernel ^Kn (^, for a suitable choice of c > 0 , and obtain asymptotics for it as n —>■ oo and 
t — 7 - 0. We will see that a transition in the asymptotics for the kernel takes place if t is of the 
order n“^, or equivalently if the distance between the two singularities in the weight is of the 
order n~^, which is of the same order as the typical distance between consecutive eigenvalues. 

Local scaling limits near the origin for the correlation kernel Kn are well understood for t ^ 0 
fixed and for t = 0, as n —?■ oo. If t 7 ^ 0 is independent of n, then, as a slight generalization of 
results obtained in [ 11 , 12 ], we have for u, u G M that 


lim 

n^oo 


, , , K, 
V’y(0)n 


V’y(0)n’ '0y(O)n 




sin7r(ri — v) 
tt{u — v) 


( 1 . 21 ) 


The convergence is uniform for u, v in compact subsets of M. If f = 0, it was proved in [27] that 
for u, V in bounded subsets of (0, 00 ) that 


lim -—— 
n^oo 1 py[ 0 )n 


Kn 


'ipv{0)n' il!vi0)n 


= K 


Bessel 


{u,v) 


= TT\/UV 


T„+i(7ru)4_i(7ru) - J„_i(7rn)4^i(7ru) 


2{u — v) 


( 1 . 22 ) 


where is the Bessel function of order v (see [29] for a reference on Bessel functions and other 
special functions which appear throughout the paper). The convergence is uniform for u, v in 
compact subsets of ( 0 ,oo). 

We obtain large n asymptotics for the correlation kernel near the origin as n —)• 00 for t small. 
Of particular interest will be the double scaling limit where n —)• 00 and simultaneously t —?■ 0 in 
such a way that nPt tends to a non-zero constant. This will lead us to a new family of limiting 
kernels which are built out of functions associated to the Painleve V equation. 


Theorem 1.7. Let V be real analytic, satisfying (1.2) and one-cut regular with 0 in the interior 
of the support [a,b] of the equilibrium measure pv We denote ify for the density of pv, given 
by (1.9), and define 

Tn,t = 167r^V’V'(0)^n^t, (1-23) 


for t real. 

1 . n —)• 00 and simultaneously t ^ 0 in such a way that Tn,t 


± 00 , we have 


lim ■ 


1 


-Kn 


u 


'0y(O)n \V'y(0)^’'0y(O)^ 

for any n, u G M, with K®™ as in ( 1 . 21 ). 


= K““(u,u), 


(1.24) 
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2. n —)• cx) and simultaneously i —)■ 0 in such a way that Tn,t 0, we have 


lim ■ 


1 


-K„ 


u 


'0\/(O)n \V'y(0)^’'0y(O)n 


= K: 


Bessel 


{u,v) 


(1.25) 


for any u,v £ M\{0}, with as in (1.22). 

3. n —)■ cx) and simultaneously t —>• 0 in such a way that Tn^t t t), there exists 
a limiting kernel such that 


lim ■ 


1 


-Kr, 


u 


'0\/(O)n \V'y(0)^’'0y(O)n 


= ]KP^(n,n;r), 


(1.26) 


for any n, n G M i/r < 0, and for any u,v G M\{±-y/r/4} i/r > 0. 


The limits are uniform for u, v in compact subsets o/M. 
Remark 1.8. The limiting kernel ^^^(n, n;r) has the form 


rPV 




<hi(7rn; r)<h2(vrn; r) — <l>i(7rn; r)<h2(vrn; r) 
2iTi{u — v) 


(1.27) 


and the functions j = 1, 2, will be defined in Section 2 for r < 0 and in Section 3 for r > 0. 
The easiest way to define <hi and $2 is via a RH problem which appeared in [4, 5]. Alternatively, 
they can be characterized as special solutions to a Lax pair which is related to the fifth Painleve 
equation. For r < 0, 4>i(n; r) and <^ 2 ( 14 ; r) are analytic as functions of n in a neighborhood of M; 
for r > 0, they are analytic functions of u in a neighborhood of M\{±^}; they are analytic 
as functions of r for r in a neighborhood of R\{0} and continuous in r at 0. For r = 0, the 
limiting kernel is a special case of a large family of limiting kernels which appeared in [2] . 

Remark 1.9. In general, the kernel is a transcendental function of u, v, r. However, for 
T < 0 and a G N, and for r > 0 and a even, the kernel is algebraic and can be constructed 
explicitly by a recursive procedure. We will present this procedure in Section 2.2 for r < 0 and 
in Section 3.2 for r > 0. 


Remark 1.10. The limiting kernel ]K^^(u, u;r) degenerates to the sine and Bessel kernels for 
large and small values of r. Indeed, we will show that 


lim ]K^^(u, u; r) = ]K®“(m, u), u,v £ 


r—>±00 


and that 


lini]KP'^(M,u;r) = K, 

T^O 


Bessel 
a 


{u, v), u,v £ M\{0}. 


(1.28) 


(1.29) 


This implies that (1.24), (1.25), and (1.26) are consistent: letting r —)• ±00 in (1.26), we 
recover (1.24), and letting r —)> 0, we get (1.25). 


1.2 Consequences and applications 
Asymptotics for Toeplitz determinants 

In [4], large n asymptotics were obtained for Toeplitz determinants 

1 p2TT 


DM) = 
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with an emerging Fisher-Hartwig singularity, depending on parameters a, (3. In the special case 
/I = 0, the weight has the form 

/,(^) = z = 0E[O,27r) 


+ 00 

for Rea > — t > 0, and where V{z) = ^ VkZ^ is analytic on an annulus containing the 

k=—oo 

unit circle. The singularities approach the unit circle as f —0, and form a Fisher-Hartwig 
type singularity if t = 0. Asymptotics for Toeplitz determinants with weight functions of this 
form were also obtained in the context of the 2d Ising model [32] , see [9] for a review of Toeplitz 
determinants and the Ising model. Theorem 1.1 in [4] states that 


OO / 

log Dnift) = nVo + ant + fc ( I 4 

k=l ^ 


^—tk 


— a- 


Vu - a- 


^—tk 


k 


+ log 


G(1 + a) 2 


r2nt 


G(1 + 2a) 


+ a^ log 2nt + 


cr„ (x) — a 


—dx + 0 ( 1 ), 


as n —)> OO and simultaneously f —)• 0 where G is Barnes’s G-function. If /3 = 0 and a = 1,2, we 
can substitute (1.17) and (1.18), here with Sn,t = “^nt, and obtain the more explicit asymptotic 
expansions 


CXJ CXJ ■ "U -/■ 

log Dnift) = nVo + ^ kVkV-k - ^(14 + V-k)e~^^ + log + ^ + o(l), 

k=\ k=l 

for a = 1, and 


log Dnift) = nVo + ^ kVkV.k - 2 ^(Vfc + + 4f 


k=l 


k=l 


— 2 log (2sinh^ t) + log (sinh^ nt — re^f^) + o(l). 


for a = 2, as n —)• OO and t —?■ 0. With more effort, we can obtain explicit asymptotic expansions 
for any integer a. 

In [5], asymptotics were obtained for Toeplitz determinants Dnift) with merging Fisher- 
Hartwig singularities, depending on parameters ai, a 2 , /3i, h- Setting ai = 02 = a/2 and 
/3i = /32 = 0, the weight has the form 

f^{z) = 

where Re a > — t E (0, vr), and with V again analytic on an annulus containing the unit circle. 
Then, the weight function has two Fisher-Hartwig type singularities if t > 0 which merge to 
a single one as t —>• 0. In [5, Theorem 1.5], the following result is stated: 

log Dnift) = logDnifo) + J - f crcK(s) “ yj ds 

oP" sin /■ o/ 

- y log — - -iVied^) + F(e-*) - 2F(1)) + o(l) 
as n —)■ OO and t —)■ 0. The function acK is related to af by the formula 


o-ck(s) 


= is) 


o? as 
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For a = 2, we can substitute (1.18) (with Sn,t = —2int) and obtain 

OO 

log Dn{ft) = nVo + kVkV-k — 2 log 2t sin t + log (n^t^ — sin^ nt) 

k=l 

- {V (e'') +F(e-'') -2^o) +o(l), 

as n —>• OO and t —?■ 0. With more effort, we can obtain explicit asymptotic expansions for any 
even a. 

Extreme values of GUE characteristic polynomials 

Let Lf be a random n x n GUE matrix, normalized such that the joint probability distribution 
of the eigenvalues is given by 

n 

n j=^ 

where = Zn{a = 0,V (x) = 2x^) is the normalizing constant. The limiting mean eigenvalue 

density is then supported on [—1,1] as n —)• oo. Define the (random) characteristic polynomial 

Pn{x) = det(x/ — H). 


The average of products of powers of characteristic polynomials of the form E n \Pn{uj)\°' , 

Vj=i ) 

given mi,M 2 £ can be expressed as 

^\Yl / n (^1 “ n 1^1 “ - u2re-^^^Uxj. (i.so) 

\j=l J j=l 

The following asymptotic results as n —>■ oo were obtained by Krasovsky [25] (see also [3, 16, 
21] for aj’s integers): 


E f n I = F{n, (ui,a,)ii) (^ + ^ ( 


logn 


n 


where k = 1,2 and 


(1.31) 


F(n,(u,,a0ii) =nC(a,)(l-n2)^'/'(n/2)“.^e(2“|-'-2'°g2)a,n H (2|^*-^.l) 

j=l 

2a2 G{a + 1)^ 


— 20140 ,^ 


C{a) = 2^ 


G{2a + iy 


and G is the Barnes G-function. If ui and U 2 approach each other, we can use Theorem 1.1 to 

/ 2 \ 

obtain asymptotics for E n \Pn{uj)\^'^^ • Indeed, by (1.30) and (1.4), it follows that 

\i=i / 


E I n 1^- 

U=1 


Zn{t = ,a,V{x) = 2{x + {ui + U 2 )/ 2 y) 


^GUE 


(1.32) 
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which is seen to hold after the simple change of variables Un = Xr, 

( ^ \ 

in the numerator, we obtain asymptotics for E n \Pn{njT]- 

Vi=i / 


. Substituting (1.15) 


This observation can be used to obtain information about extreme values of \Pn{x)\ for 
large n. This problem was investigated in [20], and in this context the authors needed large n 
asymptotics for integrals of the form (see in particular [20, Section 2]) 





n - on lim (|:Elog|Pfc(nj)|) 

e p[Uj)dUj, 


(1.33) 


with 0 G [0,1], p strictly positive and continuous on {—6,6), and a > 0. The authors note that 
one can differentiate (1.31) for k = 1 with respect to a and evaluate at a = 0 to obtain 


lim —E(log |Pn(w)|) = u — 1/2 — log2, 

n^oo Tl 


(1.34) 


since C{0) = 1, C'(0) = 0. In order to obtain asymptotics for as n —)• oo, one needs 

f 2 \ 

asymptotics for E ]/[ |P„(rij)|" , both in the region where ui and U 2 are bounded away from 

Vi=i / 

each other (in this region, we can use Krasovsky’s result (1.31)) and in the region where ui 
and U 2 are close to each other (in this region, we need to use our expansion (1.15)). 

As a corollary to Theorem 1.1 we have 


Corollary 1.11. Let 0 < 6 < 1, a > 0, and let p be a strictly positive continuous function on 
{—6,6). Then, 


In{B,p) 


C'i(a)n"^/^(1 + o(l)) fora‘s <2, 
< C' 2 nlogn(l + o(l)) for 0 ^ = 2, 
_C' 3 (a)n"^“^(l + o(l)) fora^>2, 


as n —)• oo, where 


a(l + arJ_J_e |iii-«2|»V2 


p{ui)p{u2)duidu2. 


C2 = 2 


G{1 + v/2)2 J_g 


(l — P{'^)‘^du, 


C.{a) = /" (1 - u‘) "-Pplufdu 

G{2a + l)J_g^ 


[Ml. 


-2*V+(s)-a2 


ds — iav dv. 


Remark 1.12. For o? < 2, the main contribution in the integral In{6,p) comes from the 
asymptotics (1.31) in the region of integration where ui and U 2 are bounded away from each 
other. For > 2, the main contribution comes from the asymptotics in Theorem 1.1 in the 
region of integration where \ui — tt 2 | = 0{n~^). 

Proof. First, we note that (1.31) can be extended: as n —)• oo and \ui — U 2 \ M 0 sufficiently 
slowly such that n|tti — U 2 \ M oo. 


E 



/ a a\ 

F {n,ui,U2,-,-j 


(1 + 0 ( 1 )), 


(1.35) 
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with the error term uniform for to > \ui — U 2 \ > for some to > 0 and tn ^ 0 such that 
ntn —> oo. This follows by recalling the identity 


lim 
—200 




^ as a ^ 

+ Y + Y fog 



log 


G(l + f)^G(l + 2a) 
G(1 + a)4 


(1.36) 


from [5, formula (1-26)] and applying it to (1.15), which one substitutes into (1.32). 
We now split the integral in (1.33) in two parts: the integral over 


Aa = {iui,U2): |rt2| < 0,\ui - U 2 I > 

and the integral over [—6,9]'^\A0^tnj fo^ some t„ which converges sufficiently slowly to 0 as 
re—7- 00 , slower than re“^. To compute the contribution of the integral over can 

use (1.34), (1.31), and (1.35). For the contribution of the complement of we substi¬ 

tute (1.15) in (1.32) and then, together with (1.34), in (1.33). Using (1.31) to compute the ratio 
Zn{0, a, V)(Z^^^ and changing variables t = {ui — ^2)^/4, u = (ui + U 2 )l 2 , we finally obtain 


In{e,p) 




((i-«;)(i-ui))Ws 
(2|rei - re2|)“^/2 


p{ui)p{u2)duidu2 


{l-uYl^p{uf 


(1.37) 


as re —7- 00 , where 

'sn,t = —47rmVtY(0) = —8inVt\/1 — 


in terms of the integration variables t, u. Here we need that the error terms found in Theorem 1.1 
are uniform in —9 < u < 9 for the potential V{x) = 2{x + u)‘^. This is not mentioned elsewhere 
in this paper but can readily be seen to hold true by inspection of the proof of Theorem 1.1, if 
|0| < 1. If 0 = 1, contributions from the region where the singularities approach an edge point 
of the support of the equilibrium measure will have to be taken into account as well, but this is 
outside the scope of this paper. 

From (1.36) and the small s asymptotics of it follows that with v = 4reVl — u^Vi, 


exp( r* 


Xt 


fAny/l — u^VG 


2n\/l — V? Jo 


exp 


r— 2 ™ ^-1- 


o-+(s) - W 


ds — iav dv 


= 


ore 


4v/r^ 


for < 2, 

1 G{l + V2/2fG{l + 2V2) / - 1 , ^ ^ 2 o n 

- ■!=— -re logre + o(re log re) for a = 2, (1.38) 


G(l + \/2)4 


cn ^ + o(re ^) 
as re —7- 00 , with 


for o? > 2, 


c = 


2^^\^ 


Jo 


exp 


f—2iv 


CJ+(s) - W 


ds — iav dv. 


Formulas (1.37) and (1.38) yield the corollary, where it is readily seen from (1.36) that the 
integral appearing in the constant Go is well defined. ■ 
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1.3 Outline 

In Section 2, we recall a model RH problem, introduced in [4], associated to the fifth Painleve 
equation. We define the functions 4>i and ^2 in terms of the solution to the RH problem for 
r < 0, and describe their relation with the solution a~ to the fifth Painleve equation. For a 
integer, we construct the RH solution recursively. In Section 3, we describe a second model 
RH problem which was introduced in [5], we define and ^2 in terms of its solution for 
r > 0, and we relate them to the Painleve solution cr^. For a even, we construct the RH 
solution recursively. In Section 4, we recall the standard RH problem which characterizes the 
orthogonal polynomials pj , and we define the ^-function needed to analyze this RH problem 
asymptotically. In Section 5, we perform a Deift/Zhou steepest descent analysis to obtain large n 
asymptotics for the orthogonal polynomials in the case where t < 0 (i.e., the case where the 
singularities are complex and approach 0). The construction of a local parametrix near 0, in 
terms of the model problem defined in Section 2, is crucial here. At the end of Section 5, 
we will prove Theorems 1.7 and 1.1 in the case where t < 0. In Section 6, we perform a 
similar asymptotic analysis in the case t > 0 (i.e., the case where the singularities are real and 
approach 0). Here the local parametrix is constructed in terms of the model problem from 
Section 3. This will allow us to prove Theorems 1.7 and 1.1 for t > 0. 


2 Model RH problem for r < 0 

The functions $1 and $2 that appeared in the limiting kernel can be defined in terms of 
a RH problem. We need to distinguish between the case where r > 0 and r < 0. 

The first model RH problem is a special case of the one studied in [4, Section 1.3]. For our 
purposes, the parameter /? from [4] is set to zero. The RH problem depends on parameters 
a > —1/2 and s G C. The relevant case for us will be s > 0. 


RH problem for 


(a) T-: C\r- 


C' 


2x2 


is analytic, with T = T^ U r 2 U Tg U r 4 U [0,1], and 


Ti - 2+e^ 


-r-.— 1 Stt^ 

rj = - + e . I 


hq — X + 6 


4 1 


— X + 6 


oriented as in Fig. 1. 

(b) has continuous boundary values on r“\{0, 1}, which we denote by if the limit 

is taken from the left when oriented along the contour, and 'I'Z( 2 :) if the limit is taken 
from the right. We have the jump relations 



1 5 




for 2 ; G T 
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'if+iz) = ^_{z) forzGr 2 , 

^'+(^) = (^Jia for z G Fg , 

^'+(^) = ^'_(^) ^ j for G r 4 , 

'^liz) = for z G (0,1), 


where (T 3 = 


1 0 
0 -1 


is the third Pauli matrix. 


(c) There exist q, p, r independent of 2 ; (but depending on s, a) such that 'h (z) has the 
following behavior as z —>■ 00 : 


T {z)= [I+- 

' z \p 


+ 0{z exp . 


( 2 . 1 ) 


(d) The functions G{z) := T“(z)z“ 2 ‘^® and H{z) := T“(z)(z —are bounded for z near 0 
and 1 respectively, and 'h“ is bounded near 1/2. The branches in the definitions of G, H 
are such that Z 2 , (z — 1 ) 2 > 0 on the + side of (l,oo). 


For a > —1/2, it was shown in [4, Theorem 1.8(i)] that the RH problem for T” has a unique 
solution for all s > 0 . 

It will be useful to note that the following symmetry relation holds: 

+ = ( 2 . 2 ) 


It can indeed be verified that the left and right hand side satisfy the same, uniquely solvable, 
RH problem. 

For r < 0 and u G M, we define 




for ±u > 0 , 


s = 


(2.3) 


One verifies using the jump conditions for T” that 4>i and $2 are analytic in a neighborhood 
of the real line. In fact, the singularities z = 0 and z = 1 of T” are transformed to complex 
conjugate singularities u = ±zs/4 = ±z\/—r/4 for <I>i and <I> 2 - We have the asymptotic behavior 

(I)!/))) = (^ + 0 {«-‘))exp((m-^/= 7 / 4 )cr 3 )e’=’'‘”/ 2 .. (2.4) 

as u —)• ± 00 . The kernel ]K^^(m, u;r) is then given by (1.27) for r < 0. 


2.1 Lax pair and Painleve V 

There is a connecton between the RH problem for 4'“ and the Painleve V equation, which relies 
on the theory of isomonodromy preserving deformations [14]. We recall from [4] (in particular 
Theorem 1.8(ii) and Section 4.3 in that paper) that the RH solution 41 = 4'“ solves a system of 
linear differential equations 


^^(z;^) = H(z;s)4'(z;s) 
4's(z;s) = R(z;s)4'(z;s) 


(2.5) 

( 2 . 6 ) 
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where the matrices A and B have the form 



(2.7) 

B{z-, s) = Bo{s) + Bi{s)z. 

(2.8) 


The 2x2 matrices Aq, Ai, A^o, Bq, Bi, which are independent of z but can depend on s and a, 
can be parametrized as follows, 


^0 = 

= 


-V + a/2 

V 

uy 

uy{v 

V — 

- a)\ 

a/2 


(2.9) 

- - a/2 

V 

y 

-y{v - a)\ 

—V + a/2 j ’ 


(2.10) 

s 

Bi = 

1 

Ao + Ai 

(2.11) 

2^3’ 

-2^3, 

Bo = — -i. 

s 


Here u, v, and y are functions of s and a. We have the following relation between u, v and the 
functions q, r, p appearing in the asymptotic expansion (2.1) for as z —>■ oo, 


a 

V = — — q — srp, 


u = l + 


sp 


(1 — s)p + sp' 


( 2 . 12 ) 

(2.13) 


The compatibility condition implies that Ag — Bz + AB — BA = 0, and this implies 

the following system of equations for n, v, y: 


sug = SU + {a — 2v){u — 1)^, 


svg = V \ u - \ [V — a] 

u, 

sUs = y ( —2u + a + uv — ua + — — s] . 


(2.14) 

(2.15) 

(2.16) 


By eliminating v from the top two equations, we obtain a special case of the fifth Painleve 
equation. 


^SS — 1 „ T 


1 


r, ■ 1 / ■*“ 

ZU u — 1 S 


2 , (1 - «) /^« 


1 


1 


— U-- +-- 


u u{u + 1) 
s 2{u — 1) 


If we define a{s) as 

f + OO 


cr S = 


/ + 00 

v{i)d^ 


(2.17) 


(2.18) 


(it was shown in [4] that v decays rapidly as s —>■ +oo so that this integral converges), then a 
solves (1.10). 


2.2 Recursive construction of if ck G N 

For positive integer values of a, the RH solution 4'“ can be constructed explicitly. We will 
denote 'h” instead of T” in this section to emphasize the dependence on a. The crucial ob¬ 
servation here is that the jump matrices for 4'“ are periodic in a: they remain the same if 
we replace a by a -|- 2. This fact indicates that there is a Schlesinger type transformation [28] 
relating 'I'q, to ^a+ 2 - Even more, if we replace a by a + 1, the jump matrices are modified in 
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a simple way: the jumps on the diagonals Fj”, F^, Fg , FJ are the same for the functions 4'“ 
and whereas they are opposite on (0,1). Based on this observation, and on the fact 

that the solution for a = 0 is simple and explicit, we can recursively construct a solution for any 
integer a. 

For positive integer a and s G (0, +oo), define as follows. 


Xa{z\s) 


4'“(2:;s) 


(; ‘-,"j 

1 0 

(- 1 )“ 1 


for z E II, IV, 
for z E I, 

for z E III, 


(2.19) 


with regions I, II, III, IV as in Fig. I. Then, by the RH conditions for X^ satisfies the 
following RH problem. 


RH problem for X^ 

(a) Xa : C\[0,1] —)• is analytic. 

(b) For z E (0,1), we have the jump relation Xa^+{z) = (— 1)“Vq^_(z). 

(c) As z —)> oo, 

Xaiz) = (I + 0{z~^)) exp (-ycra) ■ 

(d) The functions Ga and Ha defined by 

Gaiz) = Xaiz) ( 2 - 20 ) 

Ha{z) = Xa{z) - l)^"^ (2.21) 

are analytic at z = 0 and at z = 1 respectively. The branches in the definitions of Ga, Ha 
are such that 2 ;^, (z — 1)^ > 0 on the + side of ( 1 , 00 ). 

For 0 = 0, this RH problem is easy to solve: we have 

Xo{z) = exp (-yO- 3 ^ . 

Define 

/ z \ 

Xa{z)asXa+i{z)-^as, (2.22) 

where the square root is analytic on C\[0,1] and positive for large positive z. We will explicitly 
construct Wa in terms of Ga and Ha, and thus we will have the recursive relation 

Xa+i{z) = asWaiz)-^ Xaiz)a 3 . ( 2 . 23 ) 

From the jump relations for Xa, one verifies that Wa is meromorphic in z, with singularities 
at 0 and 1, and as z —>■ oo, we have 



Wa(z) = I + 0 


(2.24) 
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By condition (d) of the RH problem for Xq,, we obtain 

1 Eq:(^;) = f ) Ga{z)z 2'^V3 G„+i(z) V 3 for z near 0, 

/ z \ 1 

Wa{z) = (-- ) Ha{z){z - l) 2 ^^a 3 Ha+i{z) V 3 for z near 1. 


z - 1 


It follows from (2.24)-(2.26) that Wa takes the form 

jir \ T , I Qoi 

WJz) = I -\ -h 


' 1 ’ 
z z — 1 


(2.25) 

(2.26) 

(2.27) 


for matrices Pa, Qa which are independent of z. Moreover, it is easily verified that Xa exhibits 
the following symmetry, 


1 


1 


Xa\^z + -J = e 2^3aiXa \^-z + -j ai, 

which yields the relation 

-e“2°'®CTiPa(Tie2°'3 = Qa- 

Condition (2.25) gives us that 

Ga{z)~^z~ 2^^a^Waiz) is analytic in a neighborhood of 0 . 


(2.28) 


(2.29) 


(2.30) 


By substituting (2.27) into (2.30) it follows that the first row of Pa is 0, and consequently, 
by (2.30), 


z2'^^Ga{0) ^z 2‘^3a3Wa{z) is analytic in a neighborhood of 0 . 
This implies 


(2.31) 


G„,2i(0)(1 + Pa, 22 ) + G„,ii(0)P„,21 = 0 , (2.32) 

e“^Go,2l(0)Pa,21 + Ga,n{0)Pa,22 = 0. (2.33) 


By the unique solvability of the RH problem for Xa, it follows that this linear system has 
a unique solution for Pa ,22 and Pa ,21 for any positive s. 

In conclusion, since we know Xo(z), we can compute Go(0) by (2.20), and from that we can 
compute Pq by (2.32), (2.33). Substituting those in (2.29), (2.27), and (2.23), we obtain Xi(z). 
Similarly, once we know Xi(z), we can compute Gi(0), Pi, and X 2 , and so on. 

For a = 1 and a = 2, in this way, we obtain the expressions 


Xi(z) 

X 2 (z) 


1 


(e^ - 1 )y^(z - l)z 
1 


1 + (e® — l)z 


-1 


—+ z(e'^ — 1 ) 


(1 — e^(2 + s^) + e^^)(z — l)z 

'1 + e^(s — 1) — 2z + e^z(2 — 2s + s^) + z^ — e^z^(2 + s^ — e^) 
e®(—1 — s + + z(2 + s + e*(—2 + s))) 

1 — e® + e®s — z(2 + s + e®(s — 2)) 
e®(—1 — s + e®) + ze®(2 + 2s + s^ — 2e®) + z^ + z^e®(—2 — s^ + e®) 


(2.34) 

(2.35) 


e 2 
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Expanding Xa as z —)■ oo, we can obtain expressions for the Painleve V functions q(s}, r(s), 
p(s) by (2.1). In particular, for a = 0,1, 2, we obtain 


— 1 + e^(—2s + e^) 
1 — e®(2 + s^) + 


go(s) = 0, ^i(«) = Q 2 (s) = 

Using the relation (see [4, equation (4.109)]) 

= sg(s) - y, 

we can compute (Tq, ct]", and this gives (1.12)-(1.14). 

2.3 The kernel as r —>• —oo 

a 

We follow the steps of [4, Section 4.1] in the analysis as r —)■ —oo. Define 
^(z) = '^-(z) ' ' 


(2.36) 


(2.37) 


(2.38) 


where the power a/2 is taken analytic except on [0,1] and tending to 1 as 2 —)• +oo. The fraction 
cancels out the jump of T” on (0,1) and the singularities at 0 and 1, and T solves the following 
RH problem. 

RH problem for 

(a) T:C\(rrur 2 UFgUrj). 

(b) T has the jump relations 


T. 


= '^'-(^) ^_g-7ria (^)" 1 

1 0 

gTriCK ^ ^ ^ 1 


T+(z) = T_(z) 

~ ~ /l „-ma ( 2-1 j 

T+(z) = T_(z)(^^ " 


for 2 G Fj", 
for z G F^, 
for 2 G Fg , 
for z G FJ. 


(c) As z —)> oo, 

T(z) = (I + 0[z~^)) exp • 


(2.39) 


Now we will deform the jump contour in such a way that it stays away from 0, 1, and 1/2, 
to a contour as shown in Fig. 2. Define 


T(z) = < 


eW^^{z)e^p ( I fz - ^ ) (73 


for z ^ A\VJ A 2 , 


e4'^3^(z) Jj exp - ^ ) CT 3 ) forzG^, 


(2.40) 


(¥)■“ 
.0 1 


e4°'3vj/(^2;) 


exp ( I U - ^ ) <73 ) for 2 ; G A 2 , 


with Ai and A 2 given in Fig. 2. Here, (^^)^°^ is analytic except on (— 00 ,0] U [1, + 00 ). Then T 
satisfies a small norm RH problem as s —)> + 00 . 
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RH problem for 

(a) 'h : C\r —>■ is analytic, with E = Ei U r 2 given in Fig. 2. 

(b) 'h has the following jumps on E: 


$+(z) 

$+(z) 


$-(.) (J 


As s —>■ + 00 , this means that 

= ^-(z)(I + O (exp (—cs)) 

on El and r 2 , for some constant c > 0. 
(c) = I + 0(z~^) as z —>■ oo. 

By standard small norm analysis, we have 


$(z) = 1 + 0 



s —)■ + 00 , 


for z G El, 
for z G r 2 . 


(2.41) 


uniformly for z G C\r. Inverting the transformations 'h i—)• 4;' i—)• 4;' and using (2.3), we obtain 

' —Aix — s \ 


‘I>2(x;-s2) ' - ^ 


—2ix 1 
~ + 2’" 


—Aix + s 


qt - 2 


1 

-1 


+x > 0. 


By (2.41), as s —)> +oo, x G M, 


4>i(7ra; —s^) 

V—e 


+ o(.-i), 


which, by (1.27), yields the sine kernel limit as s —?• +oo, 
2 N sin7r(u —u) 


rPV 


{u,v,-s^) = 


tt{u — v) 


+ 0[s ^), 


and (1.28) is proved. In fact, it is clear from this derivation that (1-28) holds also as u, v tend 
to infinity together with s, as long as u — v converges. 
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Figure 3. The jump contour T and the jump matrices for 'I'd). 


2.4 The kernel as r —0 

Ct 

We will now study the RH problem for 4'“ in the limit where s —)• 0 and prove (1.29) 

Define 

4'd)(^;s) = et'^3^- x(z), 

where 

e for Re 2 >0, 

for Re 2 < 0. 

Then 4'd) satisfies the following RH problem. 

RH problem for 

(a) Td); C\r“ —)■ is analytic, where 

5 

r = IJfr, Fj- = {z: avgz 

i=l 

F^ = {z: argz = Svr/d}, Fg = {z: argz 

r^ = {z: argz = 77r/4}, F^ = z]R\[—iis/4,is/4]. 

The orientation is as in Fig. 3. Note that the orientation is different compared to the one 
in the jump contour for T”. 

(b) Td) has continuous boundary values on r“\{0}: 

= -1^ forzerruF^, 

4 /d)(^) = ^d)(^)(^| forzGFg-urj, 

^rd) (z) = (z)e™'"3 for zGT^. 

(c) Td) has the following behavior as z ^ oo: 

4'd)(^) = {l + 0{z-^))e^^^^xiz). 


= 7r/4}, 
= 57r/4}, 



(2.42) 


(2.43) 


(2.44) 
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Figure 4. The contour Tm and the jump matrices for M. 


(do) has the following behavior as z —?• —is/A: 

0(|z + ^|-“/2) 

VO(|z+f|"/2) 0(|2+||-«/2) 

(dl) has the following behavior as z —)• is/4: 


^(i)(2) = 




(e) is bounded at 0. 


Using (2.3) and (2.42), we obtain 


4>i(x; -s^) 
4>2(x; -s^) 


= e 4'^3\J/(1)^2;; , 


1 

-1 


(2.45) 


To study 'I'(^)( 2 ;; s) as s —7- 0, we need to recall results from [4, Sections 4.2.1 and 4.2.2]. The 
asymptotics as s —)• 0 can be described as follows: we have for | 2 ;| > <5 with 5 > 0 arbitrary that 


^3 


'h^^^(z;s) = (l + C>(|slog|s||) + M{z)e^ 2 

Here M satisfies the following RH conditions. 

RH problem for M 

(a) M : C\rM —)> is analytic, where 

Tm = 

Em,! = {z-. arg(z) = 7r/4}, Tm ,2 = {z'. arg(z) 

Em ,3 = {z-. arg(z) = 57r/4}, Tma = {z: axg{z) 

The orientation is as in Fig. 4. 

(b) M has continuous boundary values on T^ViO}: 


for TRez > 0. 


(2.46) 


3vr/4}, 

77 r/4}. 


M+{z) = M_{z) 
M+{z) = M_{z) 


1 -e"™ 

0 1 

1 

0 1 


for z e Tm,i, 
for z e Tm, 2 , 
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M+{z) = M_{z) forzeFM.s, 

M+{z) = M_{z) for2;erM,4- 

(c) M has the following behavior as z —)■ cx): 

M{z) = {l + 0{z-^))e^^^\ 

The RH problem for M has an explicit solution (which is not unique, unless one adds a condition 
to control the behaviour of M near 0) given in terms of confluent hypergeometric functions 
which is used in [4, 5]. But in the special case we are dealing with, the confluent hypergeometric 
functions degenerate to Hankel and modified Bessel functions Ha \ Ha \ la and (see [29]), 
and we use an explicit formula for M similar to the one given in [27, 31]. Writing 


= -i" ‘ 




we have 


M{z) = < 


-iH^^\{z)\ 

«+o 1 




CJlfJs 


for z G I, 
for z G II, 


-77^) (e-z) -iH\^>{e-^z) 


(e"z)5^R: 


( 1 ) 


(2.47) 


o+i 


-iz-2 


-f(a+|)(T 3 for z Gill, 


for z G IV, 


where the square roots and the Hankel and Bessel functions have branch cuts on (—oo, 0]. 
By (2.46) and (2.45), we have, for ±x > 0, 


(^_\j +C’(l«|log|s|)+ s 

Using (1.27) and substituting in (2.47), one obtains that as r —)• 0; 


1 


0 . 


rPV 


(u,u;r) = 


'Ky/UV 

8(u — v) 




- + ^ + C’(|rl^/2loglrl) +C>(lrl5+“). 

The identity Ha\z) + Ha\z) = 2Ja{z) yields our result: as r —)• 0 we have 


^PV( 


= 2 ^ T’y) {ja-^^i^^)Ja+li^u) - Ja+li^v)Ja-li^u)) 


2{u — v) 

+ 0{\T\^^Hog\T\) + 0 {\t\'2+'^). 


(2.48) 


This proves (1.29). 







Random Matrices with Merging Singularities and the Painleve V Equation 


21 



3 Model RH problem for r > 0 

The model RH problem relevant for the case r > 0 is a special case of the one found in [5] (up to 
a simple transformation), and differs from the one for r < 0, although there are similarities. For 
instance, the following RH problem, is also connected to the Painleve V equation, and an explicit 
solution can be constructed for even a. (but not for odd a, which was possible when r < 0). In 
the RH problem below, the relevant values for the parameters are a > —1/2 and s G (0, —ioo). 


RH problem for 

(a) T+ : C\r+ —is analytic, with r+ = P^ U P^ U Pj/ U Pj U [0,1], and 


I TTl I 

P^ = 1 + 




I TTl I 

T+ = l + e-^R+, 


oriented as in Fig. 5. 

(b) P"*" has continuous boundary values Pji on r^\{0,1}. We have the jump relations 


( 1 _p-ma 

0 1 

( 1 _p7r2Q:\ 

0 1 J 

'h+(z) = Pl(z) J) 

= J) 

'hi(^) = d/Kz) (J 


for z G Pj*", 
for z G P^, 
for z G Pj/, 
for z G r 4 , 
for z G (0,1) 


There exist q, r, p (depending on s, a but not on z) such that 'I'+( 2 ;) has the following 
limiting behavior as z —>■ oo: 


P+(z) 






(d) If a ^ N then the functions 

0(2) := >P+(3) (J «)2-t‘’» 


(3.1) 
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( 1 p —TTZa^X 

0 (3.2) 

are bounded for z near 0 and 1 respectively when z is in region IV, with g = 2 ism(Tra) • 
let Z 2 ^[z — 1)2 > 0 on the + side of (l,oo). If a G N then the functions 

G(z) := ^+iz) ( 3 . 3 ) 

Fr(.2) := ^r+(.2) 1)^ (^_ l)-f<^3 (3.4) 

are bounded for z near 0 and 1 respectively when z is in region IV. We let log z, log(z — 1) 
> 0 on the + side of (1, oo). 


The above RH problem is equivalent to the one from [5, Section 3]. Indeed, if we define 


^ck(C;'S) = e4‘^34,+ 



(C-^) 



(3.5) 


then Tck solves the RH problem in [5, Section 3] in the case /3i = /32 = 0, a\ = a 2 = aj^. 

In [5], it was proved that the RH problem for Tqk is solvable if a > —1/2 for every s G 
(0, —ioo), and it follows that the same is true for the RH problem for 'h'*'. 

For r > 0 and u G M, we define 




S = —ZVT 


where 


A(z) = 


e 2 \ 

Tria 1 

for Re z > 1, 

-e 2 y 



for 0 < Re z < 1 

■‘l 


rria 1 

for Re z < 0, 


—e 2 


(3.6) 


(3.7) 


and where 'I'+(z;s) has to be understood as 'I'^(z;s) for z G (0,1). Then, ^>i(u; r), <I> 2 ('u; r) are 
smooth on the real line, except at the singularities As u —)• Too, <hi and <I >2 have the 

asymptotic behavior 



similar to (2.4). 


3.1 Lax pair and Painleve V 

The results about the Lax pair and the relation with the Painleve V equation hold for 4'"^ as 
well as for T”. Indeed, T = T+ satisfies the Lax pair (2.5)-(2.6), where A and B can be 
parameterized as in (2.7)-(2.11). The functions v, u defined in (2.12), (2.13) in terms of q,r,p, 
solve equations (2.14), (2.15), and a defined as in (2.18) solves the fi-form (1.10) of Painleve V. 
It has to be noted that the functions q{s), r{s), p{s), v{s), u{s), a{s) defined here for r > 0 are 
in general different from their counterparts in Section 2 for r < 0. 
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3.2 Recursive construction for a G 2N 


Similarly to the case r < 0, we define 




for z€l, III, 
for 2 : G II, 

for 2 G IV, 


(3.9) 


with regions I, II, III, IV as in Fig. 5. Then, it is easily verified that, for even a, Xa satisfies 
the RH problem for stated in Section 2.2. It is important to note that this is not the case 
for a odd. It follows that the formulas for X 2 and q 2 given in (2.35) and (2.36) still hold for 
T > 0. By (2.37), we obtain the formula (1.14) for ■ 


3.3 The kernel Kq as t —>■ +00 

An asymptotic analysis as s = —i^/r —zoo of the RH problem for Tck, equivalent to the RH 
problem for 'h'*', has been performed in [5, Section 5]. The following result can be extracted 
from that analysis: 

+ for XG (0,1), (3.10) 

as s —>■ —zoo. Substituting (3.10) and (3.6) in (1.27), we obtain that the Painleve kernel 
converges to the sine kernel, proving (1.28) whenever |zz —z;| remains bounded in the limit where 
r —)■ + 00 : 

]K^^(zz, u; r) = ]K®“(zz, v) + O . 


3.4 The kernel as r —0 

Ct 

Recall that the parameter s in the model problem for is negative imaginary. Define 


s ) = e4°'3 4' 


X{z) 


I, 


e 2 

— Tria 


0-3 

? 

0-3 


for —< Rez < 
for Re 2 < — , 

for Re 2 > . 


Then satisfies the following RH problem: 


(3.11) 

(3.12) 


RH problem for 


(a) : C\r'’' —)• C is analytic, where T"'' is 

as in Fig. 6, 

r+ = uLir+, 

= { 2 : arg (2 - |s|/4) = 7r/4}, 

r^ = { 2 : arg (2 + s/4) = 37r/4}, 

= { 2 : arg (2 + s /4) = 57r/4}, 

r| = { 2 : arg (2 - |s|/4) = 77r/4}, 

r^ = [-1^1/4,1^1/4], 

Fe = {z- Re 2 = |s|/4}, 

T^ = { 2 : Re 2 = — s /4}. 
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e 2 


^0-3 


e 2 


0-3 





( 2 ), 






1 -1 

0 1 

1 0 

1 1 

1 - 

1 0 




for z G U f^, 
for z G r^ u r^, 
for z G r^, 


for z G r^ u r^. 


(c) 'I'^^)(z) has the following asymptotic behavior as z —>■ 00 : 

^r(2)(^) = (j + 0(z-i))e*""3^(^). (3.13) 

(d) inherits its local behavior at —and ^ from 'I''''(z) through (3.1)-(3.4) and (3.11). 
When s —)• —iO it was shown in [5, Section 6.2] that, for z bounded away from 0, 

s)x-Hz)M-\z) = 1 + e(z, s), (3.14) 

where M is as in (2.47), e(z, s) = C>(|sp) + 0{\s\‘^~^*'^) for 2a ^ Z, and e(z, s) = 0(|s| log |s|) for 
2a G Z. 

We substitute (3.14) and (3.11) into (3.6) to find that as r —>■ 0, ±x > 0, 

(_\)+£(x.s). 

By (1.27), we obtain the Bessel kernel after similar calculations as for the case r < 0, 
proving (1.29). 


4 RH problem for orthogonal polynomials 

4.1 RH problem for orthogonal polynomials and differential identity 

The correlation kernel (1.20) for the eigenvalues in the random matrix ensemble (1.1) can be 
expressed in terms of the solution of a RH problem which characterizes the orthogonal polyno¬ 
mials with respect to the weight Wn defined in (1.3). 
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Define Y by 


Y(z; n) 


(n) 


p1”’ 


(=) 


2TTiK. 


(n) 


WI \ 
Pn [X) 


X — z 


Wn{x)dx 


(n) / N 

- 4-1 / ! Wn{x)dx 


X — z 


(4.1) 


(n) (n) 

where > 0 is the leading coefficient of the normalized orthogonal polynomial pj . It is 
well-known [15] that Y can be characterized as the unique solution of the following RH problem. 

RH problem for Y 

(a) Y is analytic in C\M. 

(b) Y has the following jump relations on M, except at the singularities ±\/t in the case where 
t > 0: 


Y+{x) = y_(x) 


1 Wn{x) 
0 1 


(c) As z —)> oo, 


(d) If* > 0 and a < 0, Y(z) has the behavior 


Yiz) 


(0(1) 0{\zTVt\'^)\ 

Vo(i) o{\zTVir))^ 


as z —)• zb-v/f. If t > 0 and a > 0, E is bounded near zt-v/f. 


We can express the eigenvalue correlation kernel Kn and the logarithmic t-derivative of the 
partition function exactly in terms of Y. 

Proposition 4.1. 

(1) Let Kn be the correlation kernel defined in (1.20). For x,y G M, we have 


Kn{x,y) 


^Wnix)Wniy) (0 (y )^+(^)(^0 


27rz 


x-y 


(4.2) 


(2) Let Zn{t) = Zn{t,a,V) be the partition function defined in (1.4). Write zq = y/t when 
t > 0 and zq = iy/^ when t < 0. For f < 0, a > — ^ and for t > 0, a > 0, we have the 
differential identity 


4logZ„(*) 



(-^o) 


y-ldY 

dz 


22 


(4.3) 
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Proof. The formula for the correlation kernel (4.2) follows after a straightforward calculation 
from (1.20) and (4.1). 

To obtain (4.3), we follow ideas from [25] and [5]. In the remaining part of this proof, we 
suppress the dependence on n in our notations and we write pj , Kj , and w instead of p ^^, , 

and Wn- We start with the well known formula 

n—\ 

Zn{t) = n 

j=0 

Using the orthogonality of the polynomials and the Christoffel-Darboux formula, it follows that 


n—1 


logZ„(t) = -2^ / Pj{x)dt{pj{x))w{x)dx 


n—1 


dt 


Kn 

J=0 1 


i=o 


/*oo () \ 

\i=0 / 


w{x)dx 


f 


dt 


l^n—l 


Kr. 


idx{Pnix))Pn-lix) - Pnix)d:^{pn-l{x))) w{x)dx 


By orthogonality, it is straightforward to check that 

— logZ„(t) = -n-^-(Ji - J 2 ), 

(it Kn—l f^n 

/ oo poo 

dtPn{x)dxPn-l{x)w{x)dx, ^2 = / da;Pn{x)dtPn-l{x)w{x)dx. 

■OO J —oo 

We start by evaluating Jp. 

/ OO poo 

Pn{x)dt{dxPn-l{x)w{x))dx + / dt{pn{x)dxPn-l{x)'w{x))dx. 

■oo J —OO 


(4.4) 

(4.5) 


(4.6) 


The rightmost term in this expression vanishes, as one can see by taking the derivative outside 
the integral and using orthogonality. By the form of the weight function w given in (1.3), we 
obtain 




4‘| = 


f 


Pn{x)dxPn-l{x)w{x) 
X=F Zo 


dx. 


(4.7) 

(4.8) 


Writing dx{pm{x))\x=y = Pm,x{y), we find that 

r(l) Pnix){Pn-l,x{^) - Pn-l,x{^Zo))w{x) , , ,, . Pn{x)w{x) 

I± = / - dx + Pn-l,x[±Zo) / - dx 

' xzfzo X =F Zo 

Pn{x)w{x) 

- dx. 

X ^ Zo 


— Pn—l,a:(iZo) 


We proceed in a similar fashion for J 2 and find that 

p„_i(x)u;(x) 


KndtKn-l , « 

J2 = -n -2-^ T^Pn 

<_1 2zo 


a 

2zo 


Pn^x 


i-zo) / 

J —( 


fC 

,a;(zo) / 

J 


Pn-l{x)w{x) 
X + Zo 


dx 


X - Zo 


dx. 


Substituting (4.7)-(4.9) into (4.4), and using (4.1), it is straightforward to obtain (4.3). 


(4.9) 
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Remark 4.2. When t > 0 and a < 0, the integrals in (4.8) are not defined. Using techniques 
similar to those in [25], one can regularize those integrals and obtain a differential identity 
similar to (4.3). We refer the reader to [5, 25] for more details. 

4.2 The gr-function and normalization of the RH problem 

We proceed with a transformation of the RH problem for Y to one which is normalized to the 
identity at infinity, and which has suitable jump conditions. This is a standard part of the 
Deift/Zhou steepest descent method [13] which we apply in order to get asymptotics for Y as 
n —)> oo with t small. Define 

T{z) = ( 410 ) 

where g is defined by 

giz) = j log( 2 ; - x)dgvix), (4.11) 

with gv the equilibrium measure given by (1.9), minimizing (1.6) and satisfying the condi¬ 
tions (1.7), (1.8). Clearly g{z) = logz -t- O {z~^) as z —)■ oo, and it follows that T satisfies the 
following RH problem. 


RH problem for T 

(a) T is analytic in 

(b) T has the following jump relation for x G M (except for x = ±y/i if t > 0), 


TYx) = T_(x) 


I.t.2 _ ^^agn{g+{x)+g-{x)-V{x)+e) 


0 


„n(g+(x)-g- (x)) 


(c) T{z) = I + 0{z ^) as z —)• OO. 

(d) If f > 0, a < 0, as z —)■ ±\/t, we have 

= o{|2=F %/*!“) 

'' Vo(i) o{|2=F %/*!“) 

If f > 0 and a > 0, T is bounded near ±-v/t. 



(4.12) 


By the fact that V is one-cut regular, we have (1.9) with h real analytic and positive on [a,b]. 
By (4.11) and the Euler-Lagrange conditions (1.7), (1.8), it can be shown that g satisfies the 
following properties which we will need later on: 


(a) g+{x)-g-{x) 

(b) g+{x)+g-{x) 

(c) g+{x)+g-{x) 


0 for X > b, 

< 27ri for x < a, 

2Tri - a){b-())^dC for x G [a, 6], 

U(x) -|- £ < 0 for X > 6 and x < a, 

< a, 

G [a,b], 
> b. 


'27rf^h(^)((^-a)(^-b)y^^d^ for x 

V(x) + £ = < 0 for X 

^-2tt C h{£,) {{i - a){i - b) f^‘^ di for x 
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Figure 7. The contour Ts. 


In part (a), the square root is positive for s £ («,&), and in part (c), the square roots are 
positive for s > b and negative for s < a, with branch cut on [a, 6 ]. The function h is analytic 
on a neighborhood U of M, and we define the function cj): U\{—oo, 6 ] —)• C as 

(f){z) = Trf h{^) - a){C - d^, (4.13) 

Jb 

with the square root analytic on C\[a, b] and positive for ^ > b, and with the path of integration 
not crossing (—oo, 6 ]. By property (a) of the function g, it follows that 

T{g+{x) - g-{x)) = 2(l)±{x) for x £ (a, b). 

Combining this with property (c) for the function g we obtain 

/2n4>+{x) I 2 _^|a\ 

r+(x) = r_(x) Q g 2 n,^_(i) J for X G [a, b]. 


By property (c) of the function g it follows that 
g+ix) + g-ix) -V{x)+i = 


—2(j){x) for X > b, 

— 24 >^{x) — 27ri for x < a. 


Combined with property (a) of g it follows that 


r+(x) = r_(x) 


\ | 3;2 _ ^\a^-2n4i{x) 

0 1 


for X < o and for b < x, 


(4.14) 


where we note that e 2 n<^+(x) _ g 2 n<^_(a:) ( 4 . 14 ). 


5 Asymptotic analysis for T as n —oo for t < 0 

In what follows, we will perform further transformations on T in order to obtain a RH problem 
for which we can compute large n asymptotics. We first complete the analysis in the case t < 0. 

5.1 Opening of the lens 

We deform the jump contour for T, which is the real line, to a lens-shaped contour around [a, b] 
as in Fig. 7. It is important that the singularities ±i\/^ lie in the region outside of the lens, and 
that the lenses lie within the region U of analyticity of h. Since we want to obtain asymptotics 
which are valid for t arbitrary small, we choose the lenses to pass through 0 , and with shape 
independent of t and n. 
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Define 


/ 


T{z) 


for z outside the lens, 



for z in the upper parts of the lens, 


for z in the lower parts of the lens. 


(5.1) 


The function {z^ —1)°‘ is chosen with branch cuts on (—zoo, —U +ioo), outside the 

lens-shaped region, and positive for real z. It is such that = \z‘^ —1\°‘ ioi z £ (a, b). 


It follows that S solves the following RH problem. 

RH problem for S 

(a) S is analytic for z G C\r 5 where = T^^top U r 5 ^bottom U M is the contour shown in Fig. 7. 

(b) S has the following jump relations on T^: 




for X £ (a, b) 



(c) 5(z) = I + O [z as z —>■ oo. 

(d) S{z) is bounded as z —>■ 0, as z —>■ a, and as z —>■ 6. 

5.2 Global parametrix for t < 0 

The jump matrices for S tend to the identity matrix exponentially fast everywhere on Tg, with 
the exception of the interval (a, b) and small neighborhoods of a, b, and 0. We construct a global 
parametrix N which solves the RH problem for S, ignoring exponentially small jumps and small 

neighborhoods of a, b, and 0. We will show later on that the global parametrix is a good 

approximation of S away from a, b, 0. 

RH problem for N 

(a) N: C\[a, 6] —)• is analytic. 

(b) N has the following jump relation on {a,b), where the orientation is taken from left to 


right: 



for X £ (a, b). 


(c) N{z) = 1 + 0 (^) as z —>■ oo. 
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Define 


N{z) = D{oo)-^No{z)D{z), 


(5.2) 


where 


No{z) = , 


1 f 


i ly \ 0 7(2;) ^J \i 1 


1 i 


-1 


7(2;) = 


z-b 


z — a 


1/4 


(5.3) 


withf j analytic on C\[a, b] and positive for z > b. D{z) depends on a and t and is given by 


1/4 


D{z) = d{z)-^^, 


d{z) = exp 



°)(g - '>))*■'" 

2tt 


( t ^os\f-t\dC \ 

\Ja - a)(6- 0y 


(5.4) 

(5.5) 


with -^(2; — a) (2; — b) analytic in C\[a, b] and positive for z > b, and with pos¬ 

itive for ^ G (djb). Then, it is straightforward to verify that N solves the RH problem for N. 
Moreover, as 2; —?> a and z ^ b, N(z) has the following behavior: 


N{z) = 0 {\z — a\ as 2; —>■ a, (5-6) 

A^( 2 :) = 0 ( 12 ; — as 2; —)■ 5. (5.7) 


To verify this, one first shows that No{z) satisfies the conditions (5.6), (5.7). It then remains 
to show that d is bounded near a and b. This can be seen by deforming the integral between a 
and 6 to a contour surrounding the interval [a, b] and applying a residue argument. 


5.3 Local parametrix at a and b 

We need to construct local parametrices near the endpoints a and b. This construction, using 
the Airy function, is standard [7, 11, 12]. The explicit formula for the local parametrices near a 
and b is not relevant to later calculations, we only need that these local parametrices exist. 

Let Ub (Ua) be an open set which is independent of n and which contains b (a). The following 
RH problems can be solved in terms of the Airy function. We do not give details about this 
construction here. 


RH problem for 

(a) : Ub\Ts —)■ is analytic on Ui,\Ts- 

(b) has exactly the same jump relations as S has for 2 ; G T^ n P;,. 

(c) P^’^\z) = (I + 0{n~^))N{z) as n —)> 00 , uniformly for 2 : G dUb and \t\ sufficiently small. 

RH problem for 

(a) P(“): Ua\Ts —)• is analytic on Ua\Ts- 

(b) P^“^ has exactly the same jump relations as S has for 2 ; G T^ n Ua- 

(c) P^°'\z) = {I + 0{n~^))N{z) as n ^ 00 , uniformly for 2 ; G dUa and |t| sufficiently small. 
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5.4 Local parametrix at the origin 

Let Uq be a sufficiently small disk around the origin, which is independent of n. We want to 
construct a local parametrix which has the same jumps as S in the disk Uq and which matches 
with the global parametrix up to order 1/n at the boundary of Uq. This matching condition has 
to be valid in the double scaling limit where n —)• oo and simultaneously t —)• 0 in such a way 
that v?t tends to a non-zero constant, or equivalently such that Tn^t —)• r < 0 , with Tn^t defined 
in (1.23). 

RH problem for P 

(a) P: Uq\Ts —> is analytic. 

(b) P has the same jump relations as S': 

PZ^{z)P+{z) = SZ^{z)S+{z) for z ^UqZ Tg. 

(c) P{z) = (/ -I- O (ra“^)) N{z) uniformly for z on the boundary BUq., as n — >• oo and simulta¬ 
neously t —)• 0 in such a way that n?t tends to a non-zero constant. 

We will construct a solution to the RH problem for P in terms of the model RH solution T” 
studied in Section 2. 


5.4.1 Construction of the parametrix 

Define A: Hq —>■ C by 


X{z) = —in X < 


(j)iz) - 




-Hz) - 


HiV^) - H-iV^) 


for Im 2 ; > 0 , 
for Im 2 ; < 0 . 


(5.8) 


By (4.13), A is a conformal map in Uq (if this one is chosen sufficiently small, but containing 
and we have 


= -X{-i^/^) = + (/>(-i\^)). 


and 


A'(0) = mr'ipviH- 


We define Sn,t such that A sends i\/^ to iSn,tH to —iSn,tH- 

Sn,t = —4iA(i\/^) = -f (/>(—= —4nRe(^(i-v/^) > 0. (5-9) 


By (4.13), we have 


Sn,i = —dvrnRe / h{s){{s — a){s 



= 4n7r^/^V’y(0) + 0(n|t|^/^) = 


(5.10) 


as t —)• 0, n —)■ oo, with Tn,t as in (1.23). 

There remains some freedom to choose the contour r 5 , which has to be of the shape shown 
in Fig. 7, but is otherwise arbitrary. We choose T^ such that the upper and lower lenses are 
mapped to the straight lines Tj",..., TJ (see Fig. 3) by A. 
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Recall the definition of from (2.42). We take P of the form 

P{z) = E{z)^^^\X{zy,Sn,t)W{z). 

Here, E has to be analytic in Uq, and W takes the form 

/ ~ i)“/ 2 g-«'<i>( 2 )^'^® a'iGi for Imz > 0, 

I ((z^ — forlmz<0. 


(5.11) 


(5.12) 


It is constructed in such a way that P has the required jump conditions on F^ n I/q. The branch 
cuts of (z^ — t)“/^ are taken on the lines {—ioo,—iy/^) and {iy/^,ioo). E is an analytic 
matrix-valued function on Uq which is constructed such that P{z) approximates N(z) on the 
boundary dUo: 


E{z) = N{z)W-\z) 


Kz) - ^ 

^A(z) + 1 ^, 


VO'S 


g-a(^)^3g^<T3^ 


(5.13) 


where we have chosen the branch cuts for z in {iy/^,ioo) and —zoo) such that 


A(z) - 

,A(^)+ % 


is„,t \ 2 °'3 


= 1 -|- C>(n for Re(z) > 0 and z G dU, 


(5.14) 


as n —>■ oo and Sn,t remains bounded. Note that N{z)W ^(z) is continuous across the real line 


and that W 


-1 m- 


AR)+^ 




has no singularities at Ei^/^. No other component of E has 


discontinuities. It follows that E is analytic on Uq. 


Proposition 5.1. Let P{z) be defined as in (5.11). Then, P satisfies the RH problem for P. 

Proof. Conditions (a) and (b) are satisfied by construction. We verify (c). For z G dUo we 
have 


P{z)N-\z) = E(z)T«(A(z); Sn,t)W{z)N-\z). 


(5.15) 


It follows from (2.44) that 
^^^\X{z);Sn,t 




Kz) - 4 


rcrs 


A(z) + 


g-iA(2)a3g^<73 =/ + c>(n“^), 


(5.16) 


as n —>■ oo, where 0{n ^) is uniform for Sn,t in compact subsets of (0,oo) and z G OUq and 
branches are chosen as in (5.14). By (5.12), we find that 


P{z)N-^{z) = E{z){l + 0{n-^))E-\z). 


(5.17) 


The proposition follows if E{z) is bounded uniformly on OUq. By (5.13) it is clear that the only 
part of E which could diverge is kF“^(z)e“*'’'(^^'^3^ which contains g’3<i>(2)-*A(2)_ 'pp with this 
we note that 


^n<f>{z)-iX{z) _ g-f 


(5.18) 


which has modulus 1 because of the symmetry (j)(z) = 4>{z). Substituting this into (5.13), one 
indeed finds that E{z) = 0(1) uniformly on the boundary OUq. ■ 
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5.5 Small norm RH problem 

Define R{z) by 


R{z) 


S{z)P{z)-^ 

for z & Uq, 

5(z)P(“)(z)-i 

for 2; G Ua, 

S{z)P^^\z)-^ 

for 2; G Ub, 

Siz)N{z)-^ 

for z G C\(Po U Pa U Ub) 


(5.19) 


It follows from the conditions of the RH problems for P, P^°'\ and N that R{z) satisfies 


RH problem for R 

(a) R is analytic on C\r/j where 

= dUa u dUb u dUo u (r5\((a, b)uUouUaU Ub))). 

(b) R has the following jump relations on E^; 

R+iz) = R.{z){l + 0{n-^)) 

as n —)> oo, uniformly for 2 ; G E^, in the double scaling limit where Tn,t —>• t < 0 . 

(c) R{z) has the following asymptotics as 2 ; —)■ 00 : 

R{z) = I + 0{z-^). 

We have that R solves a RH problem, normalized at infinity and with jump matrices that are 
uniformly close to the identity matrix. This type of RH problem is often called a small norm RH 
problem, and it implies that R itself is uniformly close to /: in the double scaling limit where 
n —)■ 00 and t —)• 0 in such a way that Tn,t —>• t < 0 , we have 

R{z) = I + 0{n-^), (5.20) 

where the error term 0{n~^) is uniform for 2 ; G C\rjj. 

Remark 5 . 2 . In fact, it can be shown that the estimate (5.20) holds not only if n —)> 00 , t —)■ 0 
in such a way that Tn,t —>■ t < 0, but also if Tn,t —^ 0 and if ^ — 00 . To see this, we first 
recall from Proposition 3.1 in [4] that (5.16) holds uniformly for Sn,t £ (0,oo) (or Tn,t < 0) as 
n —?• 00 , and it follows that the same holds for (5.17). Since E is uniformly bounded, this implies 
that condition (b) of the RH problem for R is uniform for < 0 as n —>■ 00 , and thus (5.20) 
is too. 


5.6 Proof of Theorem 1.7 for t < 0 

If Tn,t —)• T < 0, we have by (5.10) that Sn,t —>■ s G (0, +00). We have that Sn,t = (0)\/^+ 

0{n~^) as n —>■ 00. 

First, recall formula (4.2) which expresses the correlation kernel in terms of Y. Expressing Y 
in T by (4.10), we obtain 


Kn{x, y) = \x^- \y^ - 




(0 l)T^\y)T+ix) 
27ri{x — y) 
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Using the identity 5 + — |U + t’/2 = —</>+ on (a, b) and substituting S into the equation, by (5.1), 
we obtain 


Kn{x,y) 


1 / 
2m{x - y) 




S^\y)S+{x) 


/g ’T’</>+(a;)|3;2 _ ^|a/2\ 

y^ncl>+{x)\x2 _t\-a/2j ■ 


(5.21) 


Now, we use the fact that S = RP for 2 ; close to the origin. Substituting this and the 
definition (5.11) of P, we obtain for x,y G Uq: 


Kn{x,y) = 1 ) (^+0 (^(yy^^n,t)E^\y) 

X R^^{y)R+{x)E+{x)^^^\x{x)-Sn,t) • 


Up to this point, no approximations have been made. Now, we use the fact that R = I + 0{n ^) 
in the double scaling limit and the analyticity of E and R, which imply that 

E~\y)R~\y)R+{x)E+{x) = I + 0(x - y), 


as X —>■ y for n sufficiently large. We find, for fixed u, v, that 

1 f _ u _ y 

ipv{0)n "■’* nipvi^y^^ nipvi^)) 

~ 2'Ki{u — v) ('^+^) iKUn)', Sn,t)^i.\x{xn)', Sn,t) (^_-^ + 0{n ), 

as n —)> 00 , t —)• 0. By (2.3) and (2.42) we have 



and it follows from (5.10) that 

= =KPV(A(x„),A(y„);-4t)+0(n-^), (5.23) 

as n —)■ 00, t —>• 0; with given in (1.27). In the case where Tn,t —)• r < 0, 

X{xn) = TTU + O {n~^), -sl^i- = Tn,t + 0{n~y (5.24) 

which proves (1.26) by continuity of the kernel. If —)> —00, the error terms in (5.24) become 
of order C)(n|t|^/^). Nevertheless we can apply (1.28) (which holds whenever u — v has a limit 
as r —)■ —00) in (5.23) to obtain 

lim ]KP^(A(x„),A(y„);-s2 i) = ]K"“(A(x„), A(y„)) + C>(s“]), 

as n —)■ 00, Sn,t 00 . This yields (1.24), since A(x„) — A(y„) = 7r(u — v) + 0(n~^} as n — 00, 
uniformly in t < tg for some sufficiently small tg > 0. Finally, if > 0, we substitute (1.29) 
into (5.23) and obtain (1.25). This completes the proof of Theorem 1.7. 
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5.7 Proof of Theorem 1.1 for t < 0 


We start from the differential identity (4.3). Following the transformations Y i—?■ T i—)• S' i—)• R 
given in (4.10), (5.1) and (5.19), and recalling the definition (5.11) of the local parametrix, we 
find the following form for Y~^Y\ Write j = 1 for Imz > 0 and j = 2 for Imz < 0. Then 


= (^B(z) + W-\z)W'(z) + 

+ ({{^(i))-i{4,C))')(A(^))) 


2,2 


(5.25) 


where ()' in general means the derivative with respect to the main argument, ()^ means the 
derivative with respect to z , and where 

B{z) = {^^^\X))~\RE)-\z){REy{z)^^^\X). 

We will show that B(z) is bounded uniformly for s G (0,oo) and n large, but we will first 
calculate the other terms in the above expression. We have 


((W-iW')(,2) + = -|e'(z) + 


a/2 ^ a/2 


z - Zo z + Zq' 


where we used 


V 

4)'±{x) + g'^ix) = 

By condition (d) of the RH problem for T” and (2.42), it follows that 


(5.26) 




as z 


Sn,t 


z- Zo 


Zo, 


(5.27) 


and 


(('I'(‘))'‘('I''‘)) 1 )a 2 (AW) 


?TLM(g-ig') (0) _ -YL + 0(1) 

Sn,t Z + Zo 


as z —)■ —zq, 


where both 0(1) terms are uniform for large n and sufficiently small 
following lemma, which we will prove below. 

Lemma 5.3. 

(1) The following identity holds: 


(5.28) 

Now, we use the 




a 

2 ' 


2 a 

(2) B(z) is hounded uniformly os n —>■ oo for t < to for some sufficiently small to- 
Equations (5.25)-(5.28) and Lemma 5.3 can be substituted into (4.3) to obtain 


|logZ„(<.aV) = -^l 


a ( 4iA'(zo) f Sn,t , (s„,t) 


Sn,t 


+ 


a 


a 


a 


2 )^ 2zo 


+ ^ (^'(-^ 0 ) - V'(zo)) + B(zo) - B(-zo) + 0(1) 


(5.29) 


as n —>■ 00 and for sufficiently small \t\. After a straightforward calculation in which one 
uses (1.11) and (5.10), Theorem 1.1 follows upon integration. 
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Proof of Lemma 5.3. Equation (2.2) implies that 

We proceed to evaluate (G“^G') 2 , 2 ( 0 ). Substituting 

'h(C) = for C in a neighborhood of 0, 

into (2.5), (2.7) and evaluating the terms of order z~^, we find that 

9,, 

{GajG-'), ,( 0 ) = - - 1 . 

Writing 

G(z) = Go(l + Giz + 0 ( 2 ;^)), 2 ; —0, 

we also have 

(G-iG')(0) = Gi, 
and, by ( 2 . 8 ), 

l-Go = BoGo, 
os 

= G^^B.Go = -^Go V3G0. 

By (5.30), (5.31) and (5.33), we obtain 




(5.30) 


ds^ Q, ' 2 

By analyzing at the behavior of 'll as s —>■ +00 in Section 2.3, one finds that as s —)> + 00 , 

(G"^G ')2 2(0) = “ 2 2 


(5.31) 

(5.32) 

(5.33) 

(5.34) 

(5.35) 


Integrating (5.34) using (2.18) and comparing (5.35) to (1.11), one finds that 
{o-'o'),,,( 0 ) = £ + “ 

It remains to show that B{z) is bounded. It follows from (5.20) and Remark 5.2 that 
R~^{z)R'{z) = 0[n~^) as n —> 00 , t —> 0. 

We will show that E and E~^E' are bounded uniformly as n —>■ 00 and t —)• 0. For z in a 
neighborhood of zq, we write 


E{z) = Ylfjiz), 

1=1 

fi{z) = D~^{oo)No{z), 
f 2 iz) = D{z)aia 3 {z + zo)~ 2 ^^ ( X{z) + 


iSn,t 


4 
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Mz) = (z- (a - n-t-3, 

f^(z} = 

Each fj, j = 1,... ,4 is uniformly bounded at ztzo as n —>• oo and t —>■ 0. This is clear for /i 
and /s- To see this for / 2 , one can perform a contour integral to find that 



Eor / 4 , it follows from the definition of A and the fact that 4>{z) = <^(z). One can do the same 
for z in a neighborhood of —zq. Thus 

B{±zo) = lim (t(^))“^(A(z); Sn,t)0(l)^^^^(A(2:); 

Z^±Z0 ^ ^ 

which is bounded as long as s lies in compact subsets of (0,+oo). To see that -B(±zo) is 
bounded as well as s —>■ oo, we substitute in T from (2.40) and use the fact that it converges 
to the identity. To see that B{±zo) is bounded as s —>■ 0 we refer to the small x behavior of T 
in [4] (formulas (4.23), (4.53), (4.61), (4.64), (4.67)). ■ 


6 Asymptotic analysis of T as n —oo for t > 0 

In this section, we will analyse asymptotically the RH problem for T in the case where t > 0, 
which means that the singularities ±\/t are real and approach the origin as t —?■ 0. Many of 
the notations here will be the same as the ones in Section 5 to emphasize the parallels in the 
analysis. We will refer to many arguments given in the previous section, so it is recommended 
to read Section 5 before this one. 

6.1 Opening of the lens for t > 0 

We open the lens as in Fig. 8. We define S in the same way as for t < 0 in (5.1), except that 
we use the contour pictured in Fig. 8 instead of the one in Fig. 7, and except for the fact that 
we define (z^ — t)" such that it is positive on (a, —\/i) and (\/i, b): 

^^2 _ _ 1^2 _ ^|Ogioarg(2-v^)gioarg(z+v7) 

with the conventions —37rj2 < arg(z — y/i) < 7r/2 and —tt/2 < arg(z + y/i) < 37r/2. 

The RH problem for S is the following. 
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RH problem for S 

(a) S is analytic for z G C\r 5 where Fs = Fs^op U r5bottom U M is the union of the lens and 
the real line, see Fig. 8. 

(b) S has the following jump relations on F^, 


S'+(x) = S-{x) 


I |3;2 _ ^|Og-2n(/)(a;) 


0 


1 


S+{x) = S-{x)\ 1 „ 

V y 

f 1 o\ 

S^{z) = S-{x) g 2 n,^(z) 




1 


S+{x) = S-{x) 


^2n4>+{x) 1^2 _ 

0 g 2 n<^_(a;) 


for X < a and b < x, 
for X G (o, —y/t) U {Vt, b), 
for z G u r5j)Q^tQjg,, 
for X G [—y/t, y/t]. 


(c) S{z) = I + O [z as 2 —>■ oo. 

(d) For a < 0, 5 has the following asymptotics as z —)• ±y/t: 

= ^ 0 ( 1 ) o{\z^vin 
VO(l) 0{\z^yft\‘^) 


s(z) = (^^^^ o(in 

^ ^ V^(i) o{i)j 


For q : > 0 the following asymptotics hold: 

as z —)■ ±y/i from outside the lense, 
S{z) = ^ ^|““) 0(1)) as z —>■ ±y/t from inside the lense. 

(e) S{z) is bounded as z —)■ a and as z —>■ 6. 


6.2 Global parametrix 

As for t < 0, we now ignore small neighborhoods of a, b, 0 and exponentially small jumps. 

The solution to the RH problem obtained in this way is constructed as for t < 0: N is given 
by (5.2)-(5.5), where the parameter t is positive in (5.5). 

The function N now has singularities at ±\/t. Applying a contour deformation argument on 
the integral in (5.5), we obtain the following behavior, 

N{z) = 0{\z ± , z —)• =F\/t. (6-1) 

6.3 Local parametrices at endpoints a and b 

The local parametrices at the endpoints a and b remain unchanged with respect to Section 5.3. 
They satisfy the same jumps as S near a and b, and as n —>■ oo, they match with the global 
parametrix N on the boundary of fixed disks around a and b. 


6.4 Local parametrix at the origin 

Let Uq be a sufficiently small disk around the origin which is independent of n. We will construct 
a function P which satisfies the following conditions. 
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RH problem for P 

(a) P: Uo\Ts —)• is analytic. 

(b) P has same jump relations as S: 

Pl^{z)P+{z) = SZ^{z)S+{z) for z G C/q n Eg. 


(c) As n —)• oo and simultaneously t —)• 0 in such a way that n^t tends to a non-zero constant, 
we have 

P{z) = {l + 0{n-^))N{z), 
uniformly for z on the boundary dUo- 


These RH conditions are exactly the same as in the case t < 0, but it has to be noted that 
the function S and its jump contour are different here than in the case f < 0. Therefore the 
construction of the local parametrix differs from the construction done before. 

In analogy with the case t < 0, define 


X{z) = —in X 


(p+iVi) + (l^+i-Vi) 

nz) - - - 

0(^) + 4>+{-Vi) 


for Imz > 0, 
for Imz < 0. 


( 6 . 2 ) 


By (4.13), A is a conformal map in Uq, 


X{Vt) = -X{-Vt) = - 4)+{-Vt)), 

and 


A'(0) = mripviO)- 

We have that A(z) is a conformal map which sends ^/i to \sn,t\/^ and —y/i to —|sn_t|/4, with 
Sn,t = -4iA(\/t) = -2n((/)+(\/t) - 4>+[-y/t)). (6.3) 

By (4.13), we have 

rVi 

Sn,t = —27rin / ipv{s)ds = —Aimiy/t'ipvi^) + = —iy/Tn,t + (6.4) 

J-Pi 

as f —>• 0, n —)■ oo, similar to the case t < 0. 

Recall the definition of in (3.11). We now fix the lens by requiring that it is mapped to 
the jump contours Ti, r 2 , Ta, and r 4 of by A. We search for P in the form 


Piz) = E{z)^^^HXizy,Sn,t)W{z), 


(6.5) 


where 


W{z) = < 


'((^2-t)"/^e-”'^(^))'"V 3 (Ti for Imz > 0, ReA(^) ^ 

((z2 - ^ ReA(z) ^ 

((z2 - for Imz > 0, ReA(z) G [-^, 

, ((z2 - for Imz < 0, ReA(z) G 
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The above construction was done in such a way that W induces the correct jump relations 
for P. E is an analytic function which has to be such that the matching condition (c) of the 
RH problem for P is valid for z on the boundary of Uq. Therefore, we recall the definition of x 
in (3.12) and let 

E{z) = 

Proposition 6.1. P{z), defined as in (6.5), satisfies the RH problem for P. 

Proof. Condition (a) and condition (b) hold by construction. We proceed to prove (c). By the 
definition of E{z), we have 

P{z)N-\z) = i?(z)'I/(2)(A(z); Sn,t)e-^^^^^’^^x{mr^E-\z). 

For z G dUo, it follows from (3.13) that 

T(2)(A(z);s„,i) = (/ + 0(n-i))x(A(z))e'^(^)'^3 

as n —)> oo, where the 0{n~^) is uniform for Sn,t in compact subsets of (0, —ioo) and X{z) suffi¬ 
ciently large. Using the same calculations as in the proof of Proposition 5.1 we find that E 
is bounded on dUo and so the result follows in the same manner as in the proof of Proposi¬ 
tion 5.1. ■ 

6.5 Small norm RH problem 

A small norm RH problem can be constructed in a similar way as in Section 5.5: define 

S{z)P{z)~^ for z G Uq, 

\ , S{z)P^'^\z)-'^ for z G Ua, 

S{z)P^^\z)-^ for .2 G Ub, 

S{z)N{z)-^ for z £C\{UoUUaUUb). 

Then R satisfies the following RH problem, similar to the case where r < 0. 

RH problem for R 

(a) R is analytic on C\r/j where 

Tr = dUa u dUb u dUo u (r5\((a, b)uUouUaU Ub)). 

(b) R has the following jump relations on F/j: 

R+{z) = R-{z){l + 0{n-^)) 

as n —)• oo, uniformly for 2 G F^, in the double scaling limit where Tn,t —>■ t > 0. 

(c) R{z) has the following asymptotics as 2 —)■ oo: 

R{z) = I + 0{z-^). 

In the double scaling limit where n —>• oo and at the same time t —?■ 0 in such a way that n^t 
tends to a non-zero constant, it follows that 

R{z) = I + 0{n-^), (6.6) 

uniformly for 2 off the jump contour for R. As for t < 0, one can extend this result to the cases 
where n —)• oo and t —)• 0 in such a way that n^t —)■ oo and in such a way that —)• 0. We do 
not give the details here, as the arguments are very similar to those in [5]: see equations (5.1), 
(5.17), (5.18), (5.25) in [5] for n^t —)• oo and equations (6.1), (6.28), (6.32), and (4.6) in [5] for 
n^t —)• 0 . 
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6.6 Proof of Theorem 1.7 

In the same way as when t < 0, we can find an expression for the correlation kernel in terms 
of S when t > 0: 


Kn{x,y) = 


1 


pn<j>+{y) \y2_f\a/2 

6iy) 


2TTi{x — y) \ |y2 _ i |«/2 

, / |t; 2 _ i|a/2 -n0+(a:) \ 


where 


= for xG 

forxG (a,-Vf)u(Vt,&). 


(6.7) 


( 6 . 8 ) 


Substituting the small norm RH solution R into the formula for the correlation kernel in (6.7), 
one obtains 


Kn{x,y) = 


1 


^n4>+{y) 




2m{x -y) \ |?/2 _ i|a /2 


/ p Tl(f)\ rp2> _ 4.|Oi/2 


From the asymptotic behavior (6.6) and the analyticity of R and y, we have 
E~\y)R+{y)R+{x)E{x) = 1 + 0(x - y), 

as, X ^ y for n sufficiently large. Substituting this into (6.9) along with the definition for W, 
one finds that 


Kn{x,y) = 


1 




Now, we can use (3.6) and (3.11) to express in terms of the functions <hi and ^ 2 - We have 
the relation 


4>i(n;r) 

4>2(n;r) 


, -|s|/4 < u < |s|/4, 

, nGM\[-|s|/4,|s|/4], 


0 


( 6 . 10 ) 


and after a similar calculation to the one in Section 5.6, using the fact that A'(0) = 7rn'0y(O), 
one finds that Theorem 1.7 holds for Tn,t £ (0,oo). 


6.7 Proof of Theorem 1.1 for t > 0 

We assume here that a > 0. The formulas (5.25), (5.26) which we obtained in the proof of 
Theorem 1.1 for t < 0, hold also in the case f > 0, if we replace by We obtain 

from (3.11) and condition d) in the RH problem for T'*' that as z —)> y/i, 

(('I'®)'‘('1'®)1)2,4Mx)) = + 0(1), 
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where the error term is uniform for large n and sufficiently small t. Likewise as z —)• —y/i, 

((9-®)-‘(wW):),yA{^)) = + ^(G-‘G'),yo) + 0(1), 

where the error term is uniform for large n and sufficiently small t. Thus, by (4.3), we obtain 


a ( 2iX'{zo) 


^logZn{t,a,V) = -^ 
at 2zo \ Sn,t 


{G-^G')UO)- 


2i\'i-zo) 

^n,t 




+ ^ + "^{V'i-zo) - V'izo)) + B{zo) - B{-zo) + 0(1) 


( 6 . 11 ) 


as n —)• oo and sufficiently small t. As in the case t < 0, we have 

(o->o');_,(o) = ^ +1. 

and we can use the large s asymptotics given in [5, equation (5.1)] to find that 

(G-'G')jyo) = = ^ + I - I . 

In the same manner as for t < 0, one can show that 

B{±zo)= lim (t(2))“\a(2 ;); Sn,i)0(l)4'(^)(A(2;); Sn,i), 

which is uniformly bounded for n large and s G (0, —ioo). To show this for s —>• —ioo, 
one uses equations (5.10) and (5.15) from [5]; while to show this for s —)• 0, one uses equa¬ 
tions (6.10), (6.19) and (6.23) from [5]. Substituting the formula for G and H into (6.11) with B 
bounded yields Theorem 1.1 for t, a > 0 after integration. 
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